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Abstract 

We apply a new functional perturbative approach to the calculation of the equal- 
time two-point correlation functions and the potential between static color charges to 
one-loop order in Coulomb gauge Yang- Mills theory. The functional approach proceeds 
through a solution of the Schrodinger equation for the vacuum wave functional to 
order and derives the equal-time correlation functions from a functional integral 
representation via new diagrammatic rules. We show that the results coincide with 
those obtained from the usual Lagrangian functional integral approach, extract the beta 
function and determine the anomalous dimensions of the equal-time gluon and ghost 
two-point functions and the static potential under the assumption of multiplicative 
renormalizability to all orders. 
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1 Introduction 



Finding an accurate (semi-)analytical description of the infrared sector of QCD is still one of 
the most important challenges of present-day quantum field theory. In this work we concen- 
trate on Yang-Mills theory, QCD without dynamical quarks, since it is in this sector where 
the peculiar properties of QCD, in particular the confining interaction between quarks, arise. 
Recently, much of the activity in this area has focused on the formulation and (approximate) 
solution of Yang-Mills theory in the Coulomb gauge [1J-[7J, the primary reason being that the 
Coulomb gauge Hamiltonian explicitly contains the color-Coulomb potential which furnishes 
the dominant nonperturbative contribution to the static or heavy quark potential. 

Semi-analytical functional approaches to the calculation of gluon and ghost propagators 
in the infrared, mostly using Dyson-Schwinger equations, have been successful in Landau 
gauge Yang-Mills theory [8j. In the so-called ghost dominance approximation, even very 
simple analytical solutions exist in the far infrared [U [9] . Although the consistency of these 
solutions is still under discussion, it is natural to inquire whether a similar approach could 
be useful in the Coulomb gauge. However, the breaking of Lorentz covariance through the 
Coulomb gauge condition makes the usual Lagrangian functional integral approach quite 
cumbersome in this gauge, see, e.g., Ref. [10]. For this reason, semi-analytical approaches in 
Coulomb gauge have mostly used a Hamiltonian formulation. A set of equations similar to 
Dyson-Schwinger equations is obtained from a variational principle using a Gaussian type 
of ansatzes for the vacuum wave functional in the Schrodinger representation [H |3]. In the 
ghost dominance approximation, furthermore, simple analytical solutions are available for 
the far infrared |2l yy . 

However, the status of the semi-analytical and analytical solutions in the Coulomb gauge 
is not yet entirely clear, for two reasons: first, two different solutions with an infrared scaling 
behavior (differing in the infrared exponents) have been found in both the analytical and the 
semi-analytical approaches [2]-[l], and there is as yet no theoretical guidance to what the 
physical solution should be; second, the inclusion of the Coulomb form factor (the form factor 
for the color-Coulomb potential, which measures the deviation of the Coulomb potential from 
a factorization in terms of ghost propagators) in the set of equations of Dyson-Schwinger 
type results problematic. In Ref. [3], the equation for the Coulomb form factor has been 
considered subleading compared to the equations for the gluon and ghost propagators and 
therefore treated in the tree-level approximation, while in Ref. [Hj all equations have been 
considered to be of the same order and therefore treated on an equal footing, with the result 
that solutions with infrared scaling behavior cease to exist. It should be emphasized that 
only solutions with scaling behavior can give rise to a linearly rising Coulomb potential, and 
that the latest lattice calculations also show a scaling behavior for the equal-time correlation 
functions in the deep infrared [5] . It is not clear at present how to improve the approximation 
used in the variational approach in order to arrive at a unique and consistent solution. 

On the other hand, an interesting relation between Landau and Coulomb gauge Yang- 
Mills theory has been pointed out in the ghost dominance approximation in Refs. [21 Hj: 
the equal-time correlation functions of the Hamiltonian approach in Coulomb gauge are the 
formal counterparts in three dimensions of the covariant correlation functions in Landau 
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gauge in four dimensions. Building on this analogy, a possible strategy seems to be to re- 
place the variational principle by a calculation of equal-time correlation functions in the 
Coulomb gauge and trying to formulate Dyson- Schwinger equations for the latter. In the 
present work, we take a first step in this direction: we set up a functional integral repre- 
sentation of the equal-time correlation functions (without taking a detour to the space-time 
correlation functions) that is the precise three-dimensional analogue of the usual functional 
integral representation of the covariant correlation functions in the Lagrangian approach to 
Landau gauge Yang-Mills theory. We also develop a diagrammatic representation and a set 
of Feynman rules for the equal-time correlation functions. We use this new formulation here 
in order to calculate the equal-time gluon and ghost two-point correlation functions and the 
potential for static color charges in Coulomb gauge perturbatively to one-loop order. We 
extract the one-loop beta function and determine the asymptotic ultraviolet behavior of the 
equal-time two-point functions and the static potential. We also show that our results coin- 
cide with those obtained in a Lagrangian functional integral approach |12l |13] and use the 
latter for the renormalization of the equal-time correlation functions and the static potential. 

The organization of the paper is as follows: in the next section, we determine the vacuum 
wave functional perturbatively to order g"^ from the solution of the Schrodinger equation. 
With the vacuum functional determined to the corresponding order, we turn to the calcu- 
lation of the equal-time gluon and ghost two-point correlation functions in Section 3. We 
also calculate the one-loop corrections to the static or heavy quark potential (and thus to 
the Coulomb form factor) in the same section. Although for the latter calculation we need 
to go beyond the terms that we have calculated for the vacuum functional in Section 2, 
the relevant additional contributions are quite simply determined. In Section 4, we provide 
another representation of the equal-time two-point functions by choosing equal times (zero) 
in the space-time correlation functions determined before in the Lagrangian functional inte- 
gral representation [121 [IS] of the theory. The static potential can also be obtained from a 
two-point function that arises in the Lagrangian approach. We use the alternative represen- 
tations of the two-point functions and the static potential to perform the renormalization of 
our results. We show the nonrenormalization of the ghost-gluon vertex in the same section 
and use it to determine the beta function and the asymptotic ultraviolet behavior of the 
two-point functions. We also show that the same beta function is found from considering 
the static potential. Finally, in Section 5, we briefly summarize our findings. We give some 
details on an important difference that arises between the Lagrangian and the Hamiltonian 
approach when it comes to the implementation of the Coulomb gauge in the Appendix. 

2 Perturbative vacuum functional 

It is very simple to write down a functional integral representation of the equal-time corre- 
lation functions, given that they are nothing but the vacuum expectation values of products 
of the field operators. In the Schrodinger representation of Yang-Mills theory in Coulomb 
gauge, the equal-time n-point correlation functions in (3-)momentum space have the foUow- 
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ing representation: 



(A^(pi, t = 0)4(p2, t = 0) • • ■ A/(p„, t = 0)) 

= y"D[A]5(V-A)FP(A)4«(p04(p2)---A^(p.)|^(A)r (1) 

Here, il>{A) is the true vacuum wave functional of the theory. The (absolute) square |'?/'(A)p 
then plays the role of the exponential of the negative Euclidean classical action in the 
corresponding representation of the covariant correlation functions (in Euclidean space). 
FP(A) = det[— V ■ D(A)], with the covariant derivative in the adjoint representation defined 

as 

D"''(A) =5"''V + ^^''"A^ (2) 

is the Faddeev-Popov determinant (in 3 dimensions) which forms a part of the integration 
measure for the scalar product of states in the Schrodinger representation (see Ref. |14j). 
Note that the fields Af{p) on the left-hand side of Eq. (P) are spatially transverse, p- A'^(p) = 
0. We will assume the transversality of the fields A" in all of the following formulae, which 
we could make manifest by introducing a transverse basis in momentum space. However, 
there is usually no need to do so explicitly. 

In order to write down the functional integral for the equal-time correlation functions 
explicitly, the vacuum wave functional needs to be specified. The analogy with the covariant 
theory suggests to make an exponential ansatz for this wave functional, in the spirit of the 
expansion in many-body physics We consider a full Volterra expansion of the exponent: 

,(A)^exp(-f:i/||...|| E E /SS;;S(-P..---.-P.) 

k=2 11,12, ...jik o.iA2f-,o.k 

xA-(pO---A:'^^(pfe)(27r)35(pi + ... + p.)^ . (3) 

Any normalization factor can be conveniently absorbed in the functional integration measure 
in Eq. ([T]). Terms linear in A in the exponent [k = 1) are excluded by the symmetry of the 
wave functional under global gauge transformations (in the absence of external color charges). 
Regarding notation, given that our Hamiltonian formalism is not manifestly covariant, we 
will denote the contravariant spatial components of 4-vectors by (latin) subindices. 
We insert this ansatz for the vacuum wave functional into the Schrodinger equation 

i/V(A) = Eoij{A) , (4) 

where H is the Christ-Lee Hamiltonian for Coulomb gauge Yang-Mills theory [T3j, 

^ 4 / (- PpW^fc) ^^(^^^fc) + ^^"W^^^W) 

+ ^ J d3xci3yppL_p<^(x)FP(A)(x,a|(-V-D)-i(-V^)(-V-D)-i|y,6)/(y). (5) 
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Here, 

B- = -\ ei.uFf, = (v X A'^ - I r'^A' x A^)^ (6) 
is the chromo-magnetic field, and 

p"(x)=P^(x) + r^4(x)l^ (7) 

the color charge density, including external static charges pq for later use. Note that we have 
extracted a factor g from the color charges in order to simplify the counting of orders of g 
in the rest of the paper. The notation (x, a|C|y, h) refers to the kernel of the operator C in 
an integral representation. 

In most of the following perturbative calculation, we will need the Hamiltonian only up 
to order g'^, where 



3! y (27r)3 (27r)3 {2it)^ i^Jk{pi,i - P2,i) + Ski{p2,j - Psj) + W^^k - Pi,fe)] 

X A|(pi)4(P2)^?(P3)(2vr)3(5(pi + P2 + P3) (10) 

g f d Pl d p^ Y cahefcde/j: r Si X \ \ facerbde/): r X X \ 



+ r'r\6^M - 5.fc5,7)J A^(pi)A^^(p2)A^fc(p3)Af(p4)(27r)^5(pi + P2 + P3 + P4) (H) 



The term ([9]) stems from the application of the functional derivative to the Faddeev-Popov 
determinant. In this term, stands for the number of colors, f^-'^'^f^'^'^ = Nc^""^, and 
p = p/|p| denotes a unit vector. In the absence of external charges, we get for the term ([T^ 

9^ f d^p „, ,1 „, , 

' f ''rAn-p)(N.g^ [ #^^1 (2.)^^^ (13) 



2 J (277)3 ^y y J (27r)3 (p - q)2; ^ ^ M|(-p) 



(27r)3 (2vr)3 ^ (pi + P3)2 ^ ' {pi + P^] 



X (2vr)35(pi + P2 + P3 + P4) Atip,)A%p,){2n)'——-—i2nf—j- ^ . (14) 
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In the term (fT3l) on the right-hand side, tij{q) denotes the spatially transverse projector or 
transverse Kronecker delta 

tijiq) = Sij - Qiqj . (15) 

We shall now show, explicitly up to order g"^, that there is a unique perturbative solution 
of the Schrodinger equation (jll) for the wave functional ipiA) in Eq. ([3]), if we only suppose 
that the dominant contribution to the coefficient function fk is at least of order g^'"^ for 
k >2. We will consider the case without external charges to begin with, and include charges 
Pq later on in the context of the static potential. To order g'^, the Schrodinger equation reads 

(iVc'-l) (/ (03/2(P)) i^^)'m + l I (I^A^-P) [P^- (/2(P))'] A"(P) =i?0, (16) 

where we have used that 

/2t (P, -P) = f2{p)S.,S-' = /2(-p)%5'^^ (17) 

(to be contracted with spatially transverse fields) as a consequence of the symmetry under 
the exchange of the arguments, of spatially rotational and global gauge symmetry, and of 
the fact that /2;ij(pi, P2) is only defined for pi + P2 = 0. Equation (fT6ll implies that, to the 
current order, 

/2(P) = IpI , (18) 
Eo = {Nl - 1) (1 (0J|P|) (2^)''5(0) . (19) 

Generally, the energy Eq cancels any field-independent terms multiplying the vacuum func- 
tional in the Schrodinger equation to any order in g. Eqs. (1181) and f|T9l) represent nothing 
but the well-known solution of the free [g = 0) theory. The choice of the sign in Eq. f|T8l) is 
dictated by the normalizability of the wave functional ([3]) to order (7°. As usual, (27r)^(5(0) 
is to be understood as the total volume of space. 

To the next (first) order of g, the Schrodinger equation is not much more complicated: 
it reads 



1 f d^Pi d?P2 d^P3 



3 J • „ rabc 



WT " [^jk{Pl,l - P2,l) + hl{p2,j - P3,j) + Slj{p3,k - Pl,k)] 



3! J (27r)3 (27r)3 (27r)3 

- 3|pi|/|5^,,(-pi, -P2, -p3)}A^"(pi)4(P2)^?(P3)(27r)35(pi + p2 + p3) = , (20) 

where we have already taken into account the results (1181) . fll9p and the fact that 

/al-fclPl, P2, P3) = r'73;.,.(Pl, P2, P3) , (21) 

which is a consequence of global gauge symmetry and the invariance of the vacuum wave 
functional under charge conjugation. The unique solution of Eq. (120|) with the full symmetry 
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under the exchange of the arguments is 

f^ijkiVl, P2, Ps) = - I 1 - P2,k) + 5jk{p2,i - P3,i) + Ski{P3J - Pl,j)] ■ 

iPil + IP2I + iPsI 

(22) 

This equahty, and all the following equalities with explicit spatial (Lorentz) indices, are 
proper equalities only after contracting with the corresponding number of transverse vector 
fields A, or, equivalently, after contracting every external spatial index with a transverse 
projector, for example in Eq. fl22|) the index i with tii(pi). 

Now we will move on to consider the Schrodinger equation to order g"^. On the left-hand 
side, terms with four and two powers of A appear, which have to cancel separately, and 
an A-independent term which must equal Eq to this order. We begin with the term with 
four powers of A. The quartic coupling (fTT!) in the Hamiltonian has to be cancelled by 
terms stemming from the second functional derivative in Eq. ([8]) acting on the vacuum wave 
functional, and a contribution from Eq. ([HI). To order g'^, the coefficient functions /2 of 
Eq. (|T8l) and /s of Eq. ( |22l) contribute, and the function which we will determine. As a 
result, the coefficient function in the vacuum wave functional takes the following (fully 
symmetric) form to order g'^: 

(|Pl| + ... + |P4|)/S(Pl,---,P4) 

= 9^ [f^T'^'^i^ik^ji - SuSjk) + r""" f'^'^i^ij^ki - SiiSjk) + f°''^''f'""'i5ij6ki - SikSji)] (23) 

- [fslijmiPU P2, -Pl - P2)tmn(pi + P2) fsfklniPS, P4, Pi + P2) 

+ /3;ifcm(Pl> P3, "Pl " P3)W(pi + P3)/3Sn(P2, P4, Pl + Ps) 

+ fsilmiPU Pa, -Pl - PA)tmn{pi + P4)/3;7fc„(p2, P3, Pl + P4)] (24) 



2 I ^abe ^cde X X (| Pll " |P 2|)(|P3| " |P4|) 

(P 

+ r''f'''SikSjl 



^Pl+P2)^ 



(IPiI - |P3|)(|P2| - |P4|) 



(Pl + P3)^ 

(Pl + P4)^ J 

This result for /4 is represented diagrammatically in Fig. [H Equation (123!) . divided by 
(|pi| + ... + |p4|), is interpreted as the elementary or "bare" four-gluon vertex. The role 
of the factor 2 and the signs in Fig. [T] will become clear in the next section. Equation (|2^ 
and the second diagram in Fig. [1] represent the contraction of two elementary three-gluon 
vertices, the latter being given mathematically by Eq. fl22l) . The contraction refers to spatial 
and color indices and the momenta, with opposite signs. Note that there is no "propagator" 
factor associated with the contraction (except for a transverse Kronecker delta), and there is 
a factor l/(|pi| + . . . + |p4|) for the external momenta which is unusual from a diagrammatic 
point of view. Finally, Eq. (!25|) and the last diagram in Fig. [1] describe an "elementary" 
Coulomb interaction between the external gluon lines. 
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^^'""X " +2perms.j - + 2 perms, j 

Figure 1: A diagrammatic representation of Eqs. fl23l) - fl25l) . Every diagram corresponds to 
precisely one of the Eqs. ( l23l) -(!25 l) . in the same order. The "2 perms." refer to permutations 
of the external legs. 



With this result in hand, we can go on to consider the terms quadratic in A in the 
Schrodinger equation to order g"^. The relevant contributions originate from Eqs. ([8]), (Q, 
(fT3l) . and dHl), and involve the functions /2 and f^. As a result, we obtain the following 
equation for the coefficient function /2 to order g"^: 

(/2(p)) 5 6., = [P -^\p\j (2^)3 (p_q)2 )^ 

^IIWT' ^' ^^'''^^^ " / (2l)3tprq)l '^M) ■ (26) 

The explicit result for /2 to order g"^ is 

/2IPJ IPI ^ y ^2vr)3 (p-q)2 + 2|p| 3 J (27r)3 2|p| + 2|q| ^ ' 



2|p| J (27r)3 2|p|+2|q| 



i^^jmPn ~^ ^mnQj ^njPm) ^ij(p) /'OQ^ 

^ (|p| + |q| + |p-q|)' 



N,g^ 1 /■ l + (p-q)2 (|p|-|q|) 



2 



2|p| 2 J (27r)3 2|p| +2|q| (p - q)2 



(29) 



(1 + (p ■ q) ) , (30) 



2|p| 2y (27r)3V - -/(p_q)2 
where we have used the contraction of an arbitrary tensor Tjj(p) 

^T,,-(p)t,,(p)=T*(p) (31) 

in order to extract the transverse part. We have presented the diagrams corresponding to 
Eqs. (I27])-([30D in Fig. [3 The first loop integral in Eqs. ([26]) and ([27]) results from Eq. ^ and 
is represented in Fig.[2]as a ghost loop because it stems from the Faddeev- Popov determinant. 
The following three loop diagrams in Fig. [2] are obtained by contracting two external legs in 
the diagrams of Fig. [T], see Eq. fl2^ . The last loop integral in Eq. fl2B]) . or the integral fl5U]) . 
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Figure 2: The diagrams corresponding to Eqs. fl27p - fl5Ul) . The bare propagator, the inverse 
2|p| of which appears in Eq. fl27|) . is marked with an open circle for later use. The first 
two one-loop diagrams correspond to the integrals in Eq. (l27|l . in the same order. The 
following diagrams represent Eqs. (l28l) - (!30l) . respectively. See the text for a motivation of 
the "crossed" gluon propagator notation in the last loop diagram. 

on the other hand, originates from the terms 0131) and fll4p in the Hamiltonian. Lacking 
a better notation, we distinguish this contribution from the contraction of Eq. fl25l) . the 
previous diagram, by marking the gluon propagator with a cross (because there is no term 
|q| in the denominator that would indicate the presence of an internal gluon propagator — in 
fact, the diagram may be interpreted to contain a nil-correlator, where 11 is the momentum 
conjugate to A). 

We have thus completed the determination of the (exponent of the) perturbative vacuum 
wave functional to order g"^. The result is given in Eqs. (l22l) . (I23l) - (p5l) . and (!27|) - (!30|) . to 
be substituted in Eq. ([3]). We can also extract the perturbative vacuum energy to the same 
order from the A-independent terms in the Schrodinger equation with the result 

Eo = (iVe' - 1) (1 (0j/2(p)) i^^rm (32) 

[cf. Eq. ( !T6|) ]. where Eqs. (p7|) - (!30|) have to be substituted for /2(p)- The explicit expression 
is not relevant for our purposes. We shall come back to the vacuum energy later in the 
context of the static potential in the presence of external charges. It should also be clear by 
now how to take the determination of the perturbative vacuum functional and the vacuum 
energy systematically to higher orders. 

3 Equal-time two-point correlation functions 

For the calculation of the equal-time correlation functions, we need to include the Faddeev- 
Popov determinant in the measure of the functional integral, see Eq. ([T]). For our diagram- 
matic procedure, it is very convenient to introduce ghost fields and write 



FP(A) = J D[c,c]exp j d\c"(x)[-V-D'^^(A) 



)m ■ (33) 
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In our conventions, we have exphcitly 

J d'x-c%^)[-V-B^\A)]c\^) = J ^c'^(-p)p2c'^(p) (34) 

+ ^ / ^^^^/""^i..-"''(pO-'(P2)^KP3)(2-)^5(Pi + P. + P3). (35) 

Note that pi j under the integral in Eq. ( l35l) can be replaced by —p2,j due to the transversality 
of A. ' 

We now have a representation of the equal-time correlation functions as a functional 
integral over the transverse components of A, the ghost and the antighost fields, see Eq. ([1]). 
The integration measure, which would be the exponential of the negative of the Euclidean 
action in the usual four- dimensional formulation (in Euclidean space), is now given by the 
exponential in Eq. fl33|) and the square of Eq. ([3]). Note that the vacuum functional is real 
(at least to order g"^) because the coefficient functions fulfill the reality condition 

(/^^:l(-pi, • • • , -p/^))* = /^tl(pi' • • • ' p^) ■ (36) 

We shall use the analogy of this representation with the familiar functional integral repre- 
sentation of the covariant correlation functions in the usual four- dimensional formulation for 
the perturbative determination of the equal-time correlation functions ([1]). The correspond- 
ing Feynman rules are easily identified: the (static) gluon propagator is the inverse of 2|p|, 
cf. Eq. ( IT8l) (the factor of two is due to the square of the wave functional in the measure), the 
other contributions — 2(/2(p) — |p|) and the other coefficient functions — 2/3(pi, p2, Ps) and 
— 2/4(pi, . . . , P4) determine the two-, three-, and four-gluon vertices. Furthermore, from Eq. 
we identify the free ghost propagator 1/p^ and from Eq. (153]) the ghost-gluon vertex. 

We consider the gluon equal-time two-point function (y4"(pi)y4^(p2)) (with t = in 
the arguments of the gluon fields to be understood) first. One of the contributions to be 
taken into account is the ghost loop, constructed from two ghost-gluon vertices fl35|) and 
two ghost propagators [see Eq. flM|) ]. and furthermore two static gluon propagators from 
Eq. flTH]) for the external lines. As it turns out, this contribution is exactly cancelled by 
the other contribution with the same graph "topology" which arises from contracting one of 
the two-gluon vertices, (minus twice) the first integral in Eq. (!27|) . with two external gluon 
propagators. Both contributions are represented diagrammatically in the first line of Fig. 
ini The cancellation of the ghost loop contribution from the perturbative gluon two-point 
function (to one-loop order) is interesting given that this contribution plays a major role in 
the nonperturbative approaches to the infrared behavior of the equal-time gluon two-point 
function [2]-[l]. The cancellation of ghost loops was found to be a general feature in the 
Lagrangian functional integral approach [T0l[T2]. An alternative way to see the cancellation 
in our present approach is to write the Faddeev-Popov determinant as 

FP(A) =exp [tr ln(- V-D(A))] . (37) 

The coefficient of the term quadratic in A in tr In ( — V ■ D(A)) precisely equals twice the 
first integral in Eq. (127|) and hence cancels out in the exponent. 
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Figure 3: Diagrammatic representation of the various contributions to the gluonic equal-time 
two-point function, see Eqs. (l37|) - pT]) . The propagators marked with open circles are taken 
from Eqs. (fT8|) and respectively, while the "direct" contractions without open circles 

refer to the contractions that appear in the course of the determination of the vacuum wave 
functional, see Figs. [T] and [21 so that the corresponding parts of the diagrams translate into 
(minus two times) the mathematical expressions [divided by (|pi| + . . . + |p4|)] 

and (I27j) - fl30l) . The notation for the sum of all diagrams with the same topology is 

explained in the text, following Eq. (jH]). 

Next we turn to the tadpole contribution which is obtained from the elementary four- 
gluon vertex extracted from Eq. (J23l) appropriately contracted with three static gluon prop- 
agators. Again, there is a second contribution with the same "topology" given by the two- 
gluon vertex from the last integral in Eq. (127]) contracted with two external propagators, 
cf. the second line in Fig. [31 The sum of these two contributions to the gluon equal-time 
two-point function is (with p = pi) 

_ 2N^g'^ A r d^q 1 1 
~ (2|p|)2 3 7 (27r)3 2|p| +2|q|% 



to be multiplied with 6"'Hij{p). 

The most complicated contribution to the two-point function comes from diagrams with 
the gluon loop topology (with two three-gluon vertices). There are three different diagrams 
of this type, represented in the third line of Fig. [31 the first from contracting two three- 
gluon vertices extracted from Eq. (1221) with four static gluon propagators (two internal and 
two external), the second from contracting the part of the four-gluon vertex given by Eq. 



2Ncg^ 4 r d^q 1 
(2|p|)3 3 7 (27r)3 2|p| + 2|q| 

2Ncg^ A r d^q 1 
""(2^37 (2^0^21^ 
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with one internal and two external gluon propagators, and the third by contracting 
the two-gluon vertex from Eq. (125]) with two static gluon propagators. The sum of these 
contributions is 

(2|p|)^ 







2|p| +2|p-q| 




(27r)3 (|p 


+ 


q 


+ 1 


P - q|)^2|q| 2|p- q| 



X 



{SkmPn + SranQk " SnkPm) tmr{P " q) tns{<l) {^IrPs + ^rsQl " ^slPr) tkl{p) (39) 



(again, to be multiphed with S tij{p), and p = pi). The tensor structure in this expression 
is invariant under the transformation q — p — q, a fact we can use to replace 2|p| + 2|p — q| 
in the numerator with 2|p| + |q| + |p — q|. The tensor structure itself can be simplified by 
performing the contractions explicitly. A straightforward, but somewhat tedious calculation 
gives 

{SkmPn + SmnQk - SnkPm) tmr{P - q) tns(q) {SlrPs + SrsQl " SslPr) tfc/(p) 

^(l-(p.c,n(2p- + 2q'+ P'f/_yy (40) 

At last, we turn to the contributions that involve the (non-Abelian) Coulomb potential. 
There are, again, three such terms, represented in the last line of Fig. [31 the first from the 
four-gluon vertex derived from Eq. (!25|) contracted with three static gluon propagators, and 
the other two using the two-gluon vertices corresponding to the two integrals (129|) and (!30l) 
contracted with two gluon propagators each. The sum of these terms is 



2N,g' 1 f d'q 1 + (p ■ q)^ (|p| - |q|)^ 2N,g' 1 fjN_(...^ ^.2^ lpl " lq 



(2|p|)3 2y (27r)3 2|q| (p - q)^ (2|p|)3 2 J (27r)3 ^ - ^ (p - q 

2N^g^ 1 f d^q 1 + (p • q)MpP - |qp 



2 



(2|p|)3 2 7 (27r)3 2|q| (p - q)^ 



(41) 



to be multiplied with 5 tjj(p) as before, and p = pi. On the left-hand side of Eq. fHT]) . 
we have added up the contributions from the contraction of the four-gluon vertex and the 
two-gluon vertex in Eq. (!29l) to give the first loop integral. The left-hand side of Eq. (14T|) is 
represented diagrammatically as the right-hand side of the last line in Fig. [31 

Putting it all together, the result for the equal-time gluon two-point function is, to order 
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(cc) 



Figure 4: 



Diagrammatic representation of Eq. (H5ll . 



+ 



2Ncg'^ 1 



/ 



c?^g 1 + (pi ■ q)^ Pi - 
(277)3 2|q| (pi-q)2 



5'^''t,,(pi)(27r)35(pi + p2). (44) 



Looking at the diagrams on the left-hand sides of the first three hnes in Fig. [3], it is clear that 
there is precisely one diagram among those of the same topology that is contructed exclusively 
from the elementary vertices ( l22i) . ( 123|) . and ( l35l) and the "bare" propagators taken from Eqs. 
( IT8|) and (IMI) . We will call this kind of diagram an "F-diagram". Interestingly, the sum of 
all diagrams with the same topology which we will refer to as an "E-diagram", can be 
constructed from the corresponding F-diagram by a formal operation that we call the "E- 
operator". It is denoted as "i?(-)" in Fig. [3] and will be illustrated considering the example 
of the third line in Fig. [21 starting from the mathematical expression for the corresponding 
F-diagram, one multiplies the integrand (of the integral over loop momentum) with the sum 
of all |k| where k runs over the momenta of all propagators in the diagram, and divides by 
the corresponding sum restricted to the momenta of the external propagators. Indeed, the 
result of this operation is given by Eq. ( !39l) in the (q ^ p — q)-symmetric form, see the 
remark after Eq. ( l39l) . 

The same rule applies to the second line in Fig. [3l or Eq. ( l38l) , only that the propagator 
that starts and ends at the same vertex has to be counted twice in the sum over (internal 
and external) |k|. Similarly, the E-operator can be used to sum the first two diagrams on 
the left-hand side of the last line in Fig. [31 We have to consider the Coulomb interaction as 
an elementary vertex given by fl25l) to this end, and count the gluon propagator that starts 
and ends at this vertex twice in the sum over |k| just as in the case of the other elementary 
four-gluon vertex. The same rule for the generation of the E-diagrams (given the elementary 
vertices and propagators) has been shown to hold up to two-loop order for the equal-time 
two-point function and to one-loop order for the four-point function in the context of a scalar 
0^ theory [IB]. Note that two contributions, the second diagram in the first line of Fig. [31 and 
the third diagram in the last line in the same figure corresponding to the first loop integral 
in Eq. ([27]) and to Eq. (!30|l . respectively, do not fit into this general scheme. 

The calculation of the equal-time ghost two-point function from the graphical rules is 
much simpler. One obtains directly 



corresponding to the diagrams in Fig. [H Note that one of the factors 1/Pi for the external 
ghost propagators cancels against the momentum dependence of the ghost-gluon vertices. 




<fq 1 - (pi ■ q)^ 
(27r)3(pi-q)2 2|q 



) 



5-^^ (27r)35(pi + p2) 



(45) 
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We shall close this section with a calculation of the static (heavy quark) potential, the 
energy for a configuration of static external color charges Pg(x). To this end, we introduce 
charges Pg(x) into the Hamiltonian, see Eqs. ([7]) and fll2p . Compared to the Coulomb term 
f|T3|) - f|T^ which was calculated in the absence of external charges, there are two new terms 
of order g'^: 

^ '^'^ p^(-p):^Pa(p)' (46) 



which is A-independent and hence only contributes to the vacuum energy but leaves the 
vacuum wave functional unchanged, and 

The latter term, when applied to the vacuum wave functional, generates the following (prop- 
erly symmetrized) expression to order (7^, 

-llWfSf^^ ^^^^ <(p04(PJ-45(P3)(2>r)'«(p. + p. + P3) , (48) 

which implies that the vacuum wave functional to order has to be modified in order to 
fulfill the Schrodinger equation with this new term. 

A term cancelling the expression (BHl) in the Schrodinger equation can only result from 
the second derivative term in the Hamiltonian [Eq. ([8])]. It is then simple to see that we 
have to add the expression 



^r"' f^-^ p»(pi)A^(p2)A^(p3)(27r)35(pi + p2 + ps) (49) 



_^ [ d^Pi d^P2 d^P3 rabc - i^^i ~ it^^i ^«rn 
2 J (27r)3 (27r)3 (2vr)3 ^ p? |p2| + Ipsr^^"^'^ ' 

to the negative of the exponent of the vacuum wave functional in order to satisfy the 
Schrodinger equation to order g^. This term describes the back-reaction of the vacuum 
to the presence of the external charges to order g"^. Observe that due to the presence of 
the external charges Pg(p), the coefficient function of y4^(pi)A^(p2) in the vacuum wave 
functional ceases to be of the form /2(p2)5"^ (5ij(27r)35(pi -|- pa). Furthermore, contrary to 
the terms found before, the contribution (l49l) is imaginary. The rest of the vacuum wave 
functional determined in the previous section remains without change. 

We now have to calculate the vacuum energy in the presence of the external charges. To 
order (7^, the result is the former one, Eq. fl32l) . without any contribution from the new term 
fH9|) . plus Eq. fH6l) which is the part of the energy that depends on the external charges and 
hence defines the potential to this order. Of course, this is just the well-known Coulomb 
potential of electrodynamics. We are really interested in the first quantum corrections to 
this "bare" potential, which are of order g^. 

In general, the vacuum energy is given by 



■ (50) 

A=0 
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which reduces to Eq. (132|) in the absence of external charges. A contribution of order 
can hence only originate from the terms in the vacuum wave functional that are quadratic 
in A. As long as we are only interested in the potential between static sources, we can 
concentrate on terms that contain precisely two powers of pq. We then start by identifying 
all the contributions to the Schrodinger equation of order that contain two powers of A 
and two powers of Pg. One of these contributions results from expanding the Coulomb kernel 

(x,a|(-V ■ D)-i(-V^)(-V • B)-'\y,b) (51) 

in Eq. ([5]) to second order in A for p = pq. The result is the term 



- f ^^^^ '^^P4 f race rbde Pl,iP2,j , nade rbce Pi-,jP'2,i 

4 V (27r)3 ■ ■ ■ (27r)3 ^ P? (pi + Ps)^ ^ ^ P?P^(pi + P4) 



X p^(pi)p^(p2)A^(p3)A^^(p4)(2vr)35(pi + . . . + p4) (52) 

on the left-hand side of the Schrodinger equation. 

Another contribution of the same type arises from the second functional derivative in Eq. 
([8]) acting (twice) on the term (H9l) . which gives the contribution 

g"^ f d^pi d?pi 1 



r d'pi d'p^ 1 ( Me. ^„ , „ ^ 

/ 7n~v^ " ' 7n~v^~ — 2\J J %lPl + P3j 
J (27r)3 (27r)3pfp^ ^ 

+ r'7''=^t.,(Pi + P4 
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|Pl + P4 




P3 
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P4I 
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Pi + P4 
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P4I 



X p^(pi)pS(p2)A^(p3)A^'(p4)(27r)35(pi + . . . + p4) (53) 

to the Schrodinger equation. The last contribution of the same type comes from the "mixed" 
term where the operator (HTll acts upon the expression (H9i) in the wave functional. The result 
is 



g^ f d^pi d^pi 1 



4 J (27r)3 (2vr)3p2p2 



|Pl + P3 




P3I 


|Pl + P3 
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PsI 



|Pl + P3 




P4I 


|Pl + P3 
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P4I 



+ r'^f'^uApi + P4) ( : , + 



|Pl + P4 




P3I 


|Pl + P4 
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P3I 



|Pl + P4 




P4I 


|Pl + P4 


+ 


P4I 



X p^(pi)p;(p2)A^(p3)A?(p4)(27r)35(pi + . . . + P4) , (54) 

to be included in the Schrodinger equation. It can be shown that no other contributions 
quadratic in A and in pq exist to order g'^. 

In analogy to the determination of the expression fHU]) from Eq. fHHl) . the three contri- 
butions flS2]) - flM|) to the Schrodinger equation are taken care of by including the following 
expression in the negative exponent of the vacuum wave functional [in addition to (122|) . 
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Figure 5: Diagrammatic representation of the contributions f l55|) (multiplied by 2) to the 
vacuum wave functional. 
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P^(Pi)p;(P2)4'(p3)^?(p4)(2vr)35(pi + . . . + P4) . (55) 



This result (multiplied by 2) is represented diagrammatically in Fig. El where we have denoted 
the "Coulomb propagator" 1/p^ as a double line. The first diagram (and its permutation) 
corresponds to the expression fl52l) . while the second diagram (plus its permutation) corre- 
sponds to the sum of the expressions fl53|) and fl5^ . From Eq. fl50|) . we find the contribution 
to the vacuum energy 
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N^g^ f (Pq 1 
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(56) 



This latter expression can be simplified by shifting q — q — p in the last two terms (in the 
round bracket). Together with Eq. fH6l) . we find for the part of the vacuum energy that is 
quadratic in the external static charge pg, 



p{p9>2) _ 9 
-^0 ~ o 



'^'^ Po(-p)^(p)p:(p) 



(2vr) 



3 ^1^ 



(57) 
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Figure 6: A diagrammatic interpretation of the static potential to order g"^. 

the following result for the static potential to order g"^ 
, 1 t Sq 1 - (p ■ q)^ 



Ncq"^ f d?q A , ((p-q)■q)^^ (|p-q|-|q| 



,2 



p2)2 j (27r)3 I (p - q)2 q2 / 2|p - q| 2|q| (|p - q| + |q|) ' 



(59) 



Note that to the proper color Coulomb potential (see, e.g., Ref. [IS]), only the antiscreening 
term, the integral in Eq. flSSl) . contributes, while the full static potential also contains the 
screening contribution (l59l) . The semi-analytical variational approaches [Sj [11] have only 
considered the proper color Coulomb potential so far. 

We can associate diagrams with the different contributions in Eqs. f l58|) - fl59!) in a natural 
way. The vertex that joins two Coulomb (double) lines and one gluon line corresponds to the 
same mathematical expression as the ghost-gluon vertex since both objects originate from 
the Faddeev- Popov operator (—V ■ D) (or its inverse). On the other hand, the vertex with 
two gluon lines and one Coulomb line translates to the expression 



IP2I 




IpsI 


IP2I 


+ 


IpsI 



^gfa,c ip^p. , (60) 



where the gluon lines carry the momenta p2 and p3, (spatial) Lorentz indices j and k, and 
color indices b and c [cf. Eq. (1491) ]. Note that the "elementary" Coulomb interaction fl25l) 
is different from the contraction of two such vertices with a Coulomb propagator. With 
these conventions, we can represent the static potential as in Fig. O The E-operator in Fig. 
|6] exclusively refers to the internal gluon propagators and thus amounts to multiplying the 
integrand with |p — q| + |q|. 

We hence have succeeded in calculating the equal-time gluon and ghost two-point func- 
tions and the static potential to one-loop order in our functional perturbative approach, with 
the results and The same results can be obtained from a straightforward 

application of Rayleigh-Schrodinger perturbation theory [T7| . Compared to these latter cal- 
culations, we have here developed a functional integral and diagrammatical approach that is 
potentially advantageous in higher-order perturbative calculations. We have also described 
a set of simplified diagrammatic rules (the "E-operator") for the determination of equal- 
time correlation functions that is expected to carry over to higher perturbative orders and is 
hoped to eventually lead to nonperturbative equations for the equal-time correlation func- 
tions analogous to Dyson-Schwinger equations. 
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4 Lagrangian approach and renormalization 



Naive power counting shows that the results of the preceding section, Eqs. (H2|) - (H5|) and 
(l58|l - (l59l) . are ultraviolet (UV) divergent and need to be renormalized. However, some of the 
denominators occuring in the loop integrals are of a different type from those that usually 
appear in covariant perturbation theory, and efficient techniques for the handling of these 
terms have yet to be developed. These remarks apply in particular to Eqs. (1431) and fl59l) . 

The equal-time correlation functions we have been calculating are a special or limiting 
case of the usual space-time correlation functions, whence we naturally obtain the represen- 
tation 

{At{p,,t = 0)A%p,,t = 0)) = r ^^(A»(pi,pi,4)A5(p2,P2,4)) (61) 

(with the space-time correlation functions written in Euclidean space-time). Since the reg- 
ularization and renormalization program has been developed for the space-time correlation 
functions, this representation is quite useful for our purposes. In the case of Coulomb gauge 
Yang-Mills theory, however, covariance is explicitly broken through the gauge condition, and 
the calculation of the space-time correlation functions in the usual Lagrangian functional 
integral approach represents a difficulty by itself. Techniques have been developed to over- 
come these difficulties and applied in Ref. [12] to the calculation of the two-point correlation 
functions to one-loop order. 

Before properly considering the renormalization of our results, we will verify that these 
coincide on a formal level with the expressions obtained via Eq. fl6T|) . taking for the space-time 
correlation functions the formulas derived in the Lagrangian functional integral approach in 
Ref. [12]. In our notation, 



(<(p.M?W>=|^-^-/p^^ (62) 



+ 



NcQ'^ f d^q {SkmPl,n + SmnQk " ^nfcPl.m) W(Pl - q) ins(q) 

{pi - 

N,g^ 1 f d^q hiipMci) (p? 4 - 



(63) 



(pjf 2 J (2vr)4 g2 (p^ _ 



2 



|5'^^t,,(pi)(27r)^(5(pi+p2), (64) 



with pI = pI +p1 4 in Euclidean space-time. Formal integration of pi^4 and p2,4 as in Eq. (1611) 
and of the component q^ of the loop momentum, most easily using the residue theorem, leads 
to our equal-time correlation function (H2l)-(Hil). In Eq. ( l62l) . we have included the tadpole 
diagram in order that the correspondence with the equal-time gluon two-point function (H2l) - 
dH]) be term by term. The tadpole diagram was not considered explicitly in Refs. 
because it vanishes in dimensional regularization. 
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The case of the ghost two-point function is even simpler, because it is instantaneous 
already in the Lagrangian approach: from Ref. [12], 



and the fact that the dependence on pi 4 and p2,4 is exclusively through the delta function 
for energy conservation implies that (c"(pi, ti)c''(p2, ^2)) contains the factor 6(ti — ^2)- Inte- 
grating over pi^4 and p2,4 or putting ti and t2 to zero then results in a factor 5(0). In this 
case, in order to reproduce Eq. fH5|) . we integrate either over pi^4 or over p2,4, which just 
eliminates the delta function for energy conservation [more symmetrically, one may integrate 
over (pi^4 -1-^2,4) instead]. Performing the integral over ^4 then converts Eq. fl^ to Eq. fllS]) . 

Although it is certainly not surprising that the Lagrangian functional integral approach 
of Refs. [ini [12], [13] gives the same results for the equal-time correlation functions as our 
Hamiltonian approach, it is also not trivial. Concerning the gauge fixing procedure, there 
is the following important difference between the two approaches: in the Lagrangian formu- 
lation, the Weyl gauge Aq = cannot be implemented in addition to the Coulomb gauge 
condition V ■ A = ^H] • Indeed, in the first-order Lagrangian formalism, integrating out the 
Ag-field rather than setting Aq to zero yields an expression in the exponent of the measure 
for the functional integral that resembles the Christ-Lee Hamiltonian [10] . The derivation of 
the Hamilton operator ([5]) by Christ and Lee [H], on the other hand, relies on the existence 
of a gauge transformation that makes any gauge field satisfy both the Coulomb and Weyl 
gauge conditions. We discuss the possibility of simultaneously implementing the Weyl and 
Coulomb gauges in the Hamiltonian and the Lagrangian approaches in the Appendix. 

Although not defined as an equal-time correlation function in our approach, it turns out 
that the static potential (l58|) - (!59!) is related to the space-time two-point function {Aq{pi)Aq{p2)) 
in the Lagrangian functional integral approach, as was first pointed out by Zwanziger [19j. 
The formal expression for the space-time correlation function is ^L2l 

/4a. ^ .fe. ■.\_(^,^c9\f A Pl,iPl,jtij{q) 

+ ^ 1 7^— ^T"::^^ ^.(pi-q)^.(q) U^'i2.)%p,+P2). (67) 




(p?)2y (27r)>i,4 



Integrating over either pi_4 or p2,4 and over the energy component of the loop momentum, 
we obtain the antiscreening contribution fl58|) to the static potential from Eq. fl66|) because 
the latter is already instantaneous. In order to find Eq. fl59|) starting from Eq. fl67j) . we 
have to put the respective other energy component, p2,4 oi pi 4, to zero in addition [this is 
not necessary in the cases of Eqs. fIBS]) and flBBl) . because there the result of integrating over 
one of the energy components is independent of the other]. For (y4g(pi, ti)y4Q(p2, t2))5 this 
procedure amounts to integrating over the relative time ti —t2, which is, in fact, intuitively 
quite appeahng for a non-instantaneous contribution to the potential between static sources. 
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We shall now use the representation (16T!) of the equal-time gluon two-point function and 
the corresponding representations of the ghost two-point function and the static potential for 
the renormalization of these equal-time correlation functions. To this end, we make use of 
the explicit expressions obtained for Eqs. fl62l) - fl67j) in Ref. [12j in dimensional regularization. 
Thus, by integrating the result for fl62|) - fl64|) according to Eq. (|6T|) . we obtain for the equal- 
time correlation function 

(5"4,,(pi)(27r)35(pi + p2) (68) 

in the limit e — 0, where ci = 3 — 2e is the dimension of space and fi an arbitrary mass scale. 
The value of the constant Ca is not relevant to our purposes, but being an integral over an 
explicitly known function of pf^/Pi, we have carefully checked that it is finite. 

From the explicit expressions for Eqs. fl65|) - fl67|) in dimensional regularization [12], we 
find directly 

ViPl) 

This procedure to regularize the equal-time correlation functions finds further support 
in the cases where the equal-time functions in the form fj42|) -f H5|) and fl58l) - fl59l) can be 
evaluated directly in dimensional regularization (in c? = 3 — 2e dimensions). For Eqs. (jUj) 
and fj45|) . identical results are obtained in both ways [17] [also trivially for Eq. P2|) and the 
loop integral in Eq. flSHl) which is just three times the one of Eq. (jlS])]. 

The results (1681) and fIBI?]) for the equal-time two-point correlation functions can be renor- 
malized in analogy to the procedures developed for covariant theories: we introduce renor- 
malized correlation functions (or correlation functions of the renormalized fields) 

^(A^(P1)4(P2)), 

^(c^(pi)c'(P2)) . (71) 
The simplest choice of the normalization conditions is 

(^«,(PiX,(P2))|p._. = ^5'^''t..(Pi)(2vr)=^5(pi + P2) , 

(4(p04(P2)) = ^ WSip^ + p,) , (72) 

at the renormalization scale k. With these normalization conditions and the results (168|) . 
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(70) 



(4(P1)4(P2)) = 



(169|) . we obtain 



^^w-' + SK^'V"^')' '''' 

to order gf^. 

The expression (!70|) for the static potential needs to be renormahzed, too. This is most 
naturally achieved by a renormalization of the coupling constant as was first suggested in 
Ref. [20], through 

9^Hp)|,=.., = ^. (74) 



see Eq. ( 1571) . Hence, from Eq. ([70 



N^g^ 11 /I k2 _ 
(47r)^ 3 \e /i"' 



(75) 



which implies for the corresponding beta function to one-loop order. 



«^'^^M'^)= (76) 

that 

P, = -^N,. (77) 

This is the well-known result from covariant perturbation theory (for Yang-Mills theory in 
covariant gauges), and has also been found in Ref. [T2] . 

For the rest of this section, we will pursue a more conventional way of renormalizing 
the coupling constant (which, however, leads to the same result). To this end, we consider 
the equal-time ghost-gluon three-point correlation function (c"(pi)c*(p2)v4^(p3)) to order 
(one loop). The calculation of this correlation function is performed in analogy with the 
determination of the equal-time two-point correlation functions in Section 3, using the result 
for the vacuum wave functional obtained in Section 2. In this particularly simple case (and 
to the order considered), the external "propagators" (equal-time two-point functions) can 
be factorized to define the equal-time proper three-point vertex r^*'^(pi, p2, ps) as 



(C (pi)c (P2)A,(P3)) = -J (^^(^(C (Pl)c (-P4)) 

X rf ^(P4, P5, Pe) (C^(-P5)C'(P2)) (4(-p6)A^(p3)) • (78) 
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Figure 7: The proper ghost-gluon vertex to one-loop order. 



The exphcit perturbative result is 
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It is represented diagrammatically in Fig. [71 Note that due to the transversality of the gauge, 
two powers of the external momenta can be factorized from the loop integrals [cf. Eq. fH5|) ] 
and, as a result, the integrals are UV finite. This phenomenon is well-known in another 
transverse gauge, the Landau gauge [21]. For future use, we note that by very lengthy 
algebra the tensor structure in Eq. flHU]) can be simplified as follows: 



Pl,l P2,n tlmiPl - q)tnr(P2 + q) [Skm{Pl,r " P3,r " ^r) 

- SmriPl,k - P2,k - '^qk) - SrkiP2,m " P3,m + ^m)] tj/fc (Ps) (27r)^5(pi + P2 + Ps) 

(pi - q) ■ (P2 + q) 



2 q pi,k + (pi ■ P2)gfc 



[pi - q)^ (P2 + q)^ 



[q- (pi - q)][q- (P2 + q)]pi,A 



+ [pi • (pi - q)][P2 • (P2 + q)]gfe} jiiyt(p3)(27r)35(pi + ps + ps) • 

We define the renormalized coupling constant in analogy with the covariant case as 



^1) 



rS(Pi,P2,P3) 



p2=p2=p2=K2 



= Z,(K)4/'(«:)rf (pi,p2,P3) 

= -igR{K)r'"pi,ktjk{P3)i2TTfS{pi + p2 + Pa) 



p?=pi=P3=«:^ 



^2) 



at the symmetric point. As a consequence, using Eq. ( 1731) and the UV finiteness of the loop 
integrals (|79l)-(|80l). 



9Rif^) = 9 



1 + — - - In — + C 



^3) 
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with a finite constant C given by (11/6)C = (4/3) Cc + (1/2)Ca + C^, where is obtained 
from the finite loop integrals in Eqs. (1791) - (ISUj) . 

For the beta function defined in analogy with Eq. (1761) . we obtain from Eq. 



/?o = -yiV,. (84) 

This beta function coincides with the one obtained in Eq. (I77p before with the renormalized 
coupling constant defined through the static potential. The integration of the renormaliza- 
tion group equation (I76|) gives the well-known (one-loop) result 

gli'^) = , X (85) 



[and the same for g\{K) (!75|) ]. It must be noted that for renormalization group improve- 
ments like Eq. (l85l) to be sensible we have to suppose that the three-dimensional formulation 
presented here is multiplicatively renormalizable to all orders in the same way as the usual 
formulation of a renormalizable covariant quantum field theory, which is not known at present 
(even the renormalizability of the Lagrangian functional integral approach to Coulomb gauge 
Yang-Mills theory has not yet been shown). Equation (1851) and the developments to follow 
are hence to some degree speculative, but it seemed of some interest to us to explore the 
consequences of the natural assumption of multiplicative renormalizability. 

With these qualifications, we go on to use a standard renormalization group argument to 
extract the asymptotic UV behavior of the equal-time two-point correlation functions. To 
this end, differentiate Eq. ( 17T]) with respect to using the K-independence of the "bare" two- 
point functions. It is then seen that the K-dependence of the renormalized two-point functions 
is determined by the anomalous dimensions (k^ c^ln Z^ c/f^^^^)- Evaluating the latter from 
Eq. (173|) and replacing g'^ in the results with gji{n), we obtain the desired renormalization 
group equations for the equal-time two-point functions, explicitly 

'^'^(4(pi)4(P2)) = I ^^^^ (4(pi)4(P2)) • (86) 

In these equations, we substitute from Eq. fl55]) for g%{K) and integrate. Using the normal- 
ization conditions (1721) for the determination of the integration constants, one obtains the 
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momentum dependence of the equal-time two-point functions: 



2|pi 



In 



K 



3/11 



In 



A2 



rt,,(pi)(27r)^(5(pi + p2) 



(4(P1)4(P2)) 



P? 



In 



4/11 



A2 



\ 



In 



^^QCD 



5^^' (27r)^5(pi + p2) 



(87) 



The momentum dependence of the "bare" two-point functions, obtained from Eq. (!87|) simply 
by multiplying with the corresponding wave function renormalization constants Za,c, is ob- 
viously the same. By solving the renormalization group equations for Za and Zc that involve 
the anomalous dimensions, it may be shown explicitly that the bare two-point functions are 
independent, as they must be. 
For the static potential, on the other hand, we immediately obtain from Eqs. (1741) and 
f l85|l [for gji{n)] the renormalization group improved result 



yiVcP^ln 



P^ 



^^QCD 



Note that the latter one-loop formula constitutes a very direct expression of asymptotic 
freedom. 

The result (1871) for the momentum dependence of the equal-time two-point functions has 
also been obtained in Ref. [22] from a Dyson-Schwinger equation for the equal-time ghost 
correlator, where the gauge-invariant one-loop running fl85|) of the renormalized coupling con- 
stant is used as an input. We briefly discuss that derivation here, adapted to the conventions 
of the present paper. 

The renormalized equal-time two-point functions are parameterized as 



and 



(4(pi)4(p2)) 



5''^ti,(pi)(27r)^5(pi + P2) 



Pi 



^9) 



(90) 



and normalized according to the conditions (1721) . The renormalized coupling constant is 
defined as before in Eq. (182|) . Then the Dyson-Schwinger equation for the equal-time ghost 
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two-point function reads [3] 

d-\Y>'') = Z.~N,gl{H) 



(27r)3 2uj{ci^) 



(91) 



Here we have approximated the full ghost-gluon vertex appearing in the exact equation 
by the tree-level vertex, as it is appropriate in order to obtain the (renormalization-group 
improved) one-loop expressions. 

In order to solve Eq. (1911) . we make the following, properly normalized, ansatzes for the 
two-point functions in the ultraviolet limit ^ ^qcd? 



a;(p^ 



In 



K 



A2 



\ 



In 



A2 



dip' 



In 



A2 



In 



v 



A2 



(92) 



with the exponents 7 and 6 to be determined. The integral in Eq. fl9T|) can then be calculated 
in the limit ^ ^qcd ^"^^ Dyson- Schwinger equation yields the relation [22] 



In" 



K 



A2 



In^ 



P^ 



A2 



K 



from which we infer the sum rule 



(47r)2 35 
7 + 25 = 1 



A2 



for the exponents as well as the identity 

1 4 



(47r)2 35 



A2 



In 



l-7-(5 



P^ 



A2 



(93) 



(94) 



(95) 



for the coefficients. Consistency of the latter relation with the well-known perturbative result 
(153]) yields the exponents 

3 4 
7 = — , 5 = — , (96) 
' 11 ' 11 ' ^ ' 

where we have used the sum rule flM]) again. We have thus regained the result of Eq. fl57]) . 



5 Conclusions 

In this work, we have accomplished a systematic perturbative solution of the Yang-Mills 
Schrodinger equation in Coulomb gauge for the vacuum wave functional following the 
method in many-body physics. This resulted in a novel functional integral representation 
for the calculation of equal-time correlation functions. We have derived a diagrammatical 
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representation of these functions, order by order in perturbation theory, where the vertices 
in the diagrams are determined from the perturbative calculation of the vacuum wave func- 
tional. The number of the vertices, which by themselves have a perturbative expansion, 
grows with the perturbative order. We have determined the equal-time gluon and ghost 
two-point correlation functions and the potential between static color charges to one-loop 
order in this way. 

The results coincide with those of a straightforward calculation in Rayleigh-Schrodinger 
perturbation theory [T7|, and also with the values for equal times of the two-point space- 
time correlation functions from a Lagrangian functional integral representation [I2l. We 
have emphasized that the latter coincidence is not trivial since the gauge fixing procedures 
in the Hamiltonian and the Lagrangian approach are profoundly different. We have also used 
the results of the Lagrangian approach to renormalize the equal-time two-point correlation 
functions and the static potential. 

With the help of the non-renormalization of the ghost-gluon vertex which we also show, 
or, alternatively, from the static potential, we can extract the running of the correspond- 
ingly defined renormalized coupling constant. The result for the beta function is the one 
also found in covariant and other gauges, Pq = — (ll/3)A^c to one-loop order. We have used 
standard renormalization group arguments to determine the asymptotic ultraviolet behav- 
ior of the equal-time two-point functions and the static potential under the assumption of 
multiplicative renormalizability to all orders, with the result that 

(In(p2/A2 



2|p 



i| 



-4/11 



Pi 
PI 

to one-loop order in the perturbative (asymptotically free) regime. 

It is clear from the presence of an infinite number of vertices in the functional integral 
representation of the equal-time correlation functions (to infinite perturbative order) that the 
corresponding Dyson-Schwinger equations contain an infinite number of terms, a very serious 
problem for the determination of an appropriate approximation scheme for nonperturbative 
solutions. The existence of simplified diagrammatic rules for the calculation of equal-time 
correlation functions via the E-operator, to be appropriately extended to all perturbative 
orders, seems to point toward the possibility of formulating similar nonperturbative equations 
with a finite number of terms. It would indeed be very interesting to repeat the type of 
infrared analysis applied before to Yang-Mills theory in the Landau gauge [H [H [9] and to a 
variational ansatz in the Coulomb gauge [TJ-[4] for such a set of equations. 
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A Gauge transformations in Lagrangian and Hamilto- 
nian formalisms 

While the Lagrangian approach to Yang-Mills theory offers some convenient features (such 
as manifestation of Lorentz invariance), the more cumbersome Hamiltonian approach yields 
equations of motion invariant under a larger set of gauge transformations. Prior to quan- 
tization, we discuss gauge invariance starting from the classical Lagrangian and Hamilto- 
nian functions, respectively. Note that we will employ standard covariant notation in this 
appendix; in particular, spatial subindices refer to the covariant components of the corre- 
sponding 4-vector or tensor. 

The Lagrangian function of the gauge sector. 



is invariant under gauge transformations of the gauge field A^{x) = A^(x)r", 

A^{x) ^ U{x)A^{x)U^{x) + ^U{x)d^U^{x) , (99) 

where U G SU{N) and [T^ T^] = /'^*<=T^ 

The Weyl gauge, Aq{x) = 0, can be found by choosing the time-ordered exponential 



^.W^Texp{-./'...„,x,0} 



(100) 



To remain in the Weyl gauge, the transformation f llOOp may be followed by time-independent 
transformations f/(x) only. We can therefore fix the Coulomb gauge, diAl^{x) = 0, at one 
instant of time but it is impossible to fix both gauges simultaneously for all times. 

In the Hamiltonian formalism, on the other hand, gauge transformations are generated by 
(first-class) constraints in configuration space [23] • To see that, supplement the Hamiltonian 
function 

H=]^jd'x [Ulix) + B^(x)) - j d'xAl{x)Dt{x)W,{x) (101) 
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by the constraints 

= nS(x) ^ , = Dt{x)Yl',{x) ^ (102) 

with some arbitrary Lagrange muhipher fields {A^(a;)}, 




(103) 



We defined Il1{x) = F^^^^x) and t)f{x) = d^'^di - gf'^^Al^x). The extended Hamiltonian 
He in Eq. (11031) is equivalent to the original Hamiltonian H since the constraints {(t>k{^)} 
vanish weakly (in the Dirac sense [23] )• The infinitesimal time evolution of the gauge field 
y4^(x, t) from to t = to + generated by He through the Poisson brackets, 




gives for two different sets of Lagrange multipher functions {A'^^(a;)} and {A'^^(x)} two differ- 
ent results and A"^, respectively. These differ to 0{5t) by 

- KM = stJ2 f d'y {K\y) - Af(y)) {A^(x,t),0t(i/)} (105) 

and are physically equivalent. Thus, the function 

^ = E / dhr-,{ymy) (106) 

fc=l,2 

generates infinitesimal gauge transformations in the (extended) Hamiltonian formalism with 
arbitrary functions t^{x) and (x). Computing the Poisson brackets in Eq. (11051) yields 

Al{x) ^ Al{x) + r'^ix) (107) 
A1{x) ^ A1{x) - bl\x)Tl{x) (108) 

The difference to the gauge transformations (!99|) in the Lagrangian formalism is that the 
time component and the spatial components of the gauge field transform independently. 
The two functions r"(x) and r2 (a;) allow for a larger set of gauge transformations than the 
single function U{x) in the Lagrangian formalism. The simultaneous fixing of Weyl and 
Coulomb gauges, which is impossible in the Lagrangian formalism, can be accomplished in 
the Hamiltonian formalism by appropriately choosing t^{x) and (x). See Ref. [Ml for the 
abelian case. Subsequently, the non-abelian gauge-fixed theory can be canonically quantized 
with projection on the physical Hilbert space [H], or with Dirac brackets [22] enforcing 
all constraints strongly. Both quantization prescriptions produce the Hamiltonian operator 
given by Eq. ([5]). 
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